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Abstract. We have developed a homogeneous model of physical chemistry to investi- 
gate the neutral-dominated, water-based Enceladus torus. Electrons are treated as the 
summation of two isotropic Maxwellian distributions— a thermal component and a hot 
component. The effects of electron impact, electron recombination, charge exchange, and 
photochemistry are included. The mass source is neutral H 2 0, and a rigidly-corotating 
magnctosphcrc introduces energy via pickup of freshly-ionized neutrals. A small fraction 
of energy is also input by Coulomb collisions with a small population (< 1%) of supra- 
thermal electrons. Mass and energy are lost due to radial diffusion, escaping fast neu- 
trals produced by charge exchange and recombination, and a small amount of radiative 
cooling. We explore a constrained parameter space spanned by water source rate, ion 
radial diffusion, hot-electron temperature, and hot-electron density. The key findings are: 
(1) radial transport must take longer than 12 days; (2) water is input at a rate of 100- 
180 kg s" 1 ; (3) hot electrons have energies between 100 and 250 eV; (4) neutrals dom- 
inate ions by a ratio of 40:1 and continue to dominate even when thermal electrons have 
temperatures as high as 5 eV; (5) hot electrons do not exceed 1% of the total elec- 
tron population within the torus; (6) if hot electrons alone drive the observed longitu- 
dinal variation in thermal electron density, then they also drive a significant variation 
in ion composition. 



1. Introduction 

Absorption of UV starlight during occultation of Saturn's 
moon Enceladus showed that it continuously ejects neutral 
H2O at a rate of w 150-300 kg s^ 1 from water-ice geysers 
located at its southern pole [Hansen et al, 2006]. Models 
suggest that the water and its chemical by-products form 
an extended neutral-dominated torus centered on the or- 
bit of Enceladus [Jurac and Richardson, 2005]. Similarly, 
Jupiter's volcanic moon Io emits a mixture of SO2 and 
S2 at a rate of ~ltons _1 , and chemical by-products pro- 
duce a plasma torus centered on the orbit of Io (see re- 
view by Thomas et al. [2004]). The Hubble Space Tele- 
scope (HST) observations by Shemansky et al. [1993] re- 
vealed that the neutral-to-ion ratio in the Enceladus torus 
(~ 10) is three orders of magnitude greater than in the Io 
torus (ss 10~ ). Compositional differences and the degree 
of ionization within these two systems can be attributed 
to their chemistry (Io's based on sulfur dioxide and Ence- 
ladus's based on water) as well as the fact that fresh ions 
are picked with five times more energy in the Io torus than 
in the Enceladus torus [Delamere et a/., 2007]. 

An important lesson learned from studying the physi- 
cal chemistry of the Io torus is that a small fraction of 
hot electrons (< 1%) play a critical role in determining 
composition. To model the Cassini UltraViolet Imaging 
Spectrograph (UVIS) data obtained during Cassini's E2 
flyby of Jupiter (October 2000 to March 2001), Steffi et al. 
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[2004a, b, 2006, 2008] adapted the Delamere and Bagenal 
[2003] Io torus model. Steffi et al. used their models to 
study radial, temporal, and azimuthal variation in mixing 
ratios (ion-to-electron density ratios) , thermal-electron den- 
sity, and thermal-electron temperature. Steffi et al. [2008] 
concluded that hot electrons are necessary for the Io torus 
energy budget and that two modulations of the hot-electron 
population are required to reproduce both the temporal and 
spatial variations in composition observed in the data, one 
modulating in Jupiter's System III longitude, the other in 
System IV. We anticipate that hot electrons are similarly 
important in Saturn's Enceladus torus. 

Delamere et al. [2007] developed a simplified oxygen- 
based model to compare the Enceladus and Io tori. They 
found that collisional heating by a population of hot elec- 
trons is much less important at Enceladus, contributing only 
0.5% of the energy to the torus, compared to 60% at Io. 
They also cited two major reasons for the discrepancy in the 
neutral-to-ion ratio between the two systems. First, newly 
created ions are picked up in the Io torus by Jupiter's mag- 
netosphere at roughly five times the energies as are those 
in the Enceladus torus by Saturn's magnetosphere. The 
higher-energy pickup ions in the Io torus warm the thermal 
electrons, which then reduces the neutral-to-ion ratio via 
impact ionization. Second, because of the high abundance 
of molecular ions compared with atomic oxygen ions in the 
Enceladus torus (e.g., Sittler et al. [2005]), Delamere et al. 
[2007] expected that molecular dissociative recombination 
(not included in their model) will therefore drive the ratio 
even higher in the Enceladus torus. From the conclusion of 
Delamere et al. [2007]: "The addition of the full water-group 
molecular chemistry will introduce an additional plasma sink 
through dissociative recombination of the molecular ions. 
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Therefore, our simplified O-based chemistry likely represents 
a lower limit for the neutral/ion ratio." To investigate the 
consequences of a water-based Enceladus torus dominated 
by molecular chemistry, we have improved on the Delamere 
model by including a comprehensive set of water-based re- 
actions and species to more accurately estimate steady-state 
densities and temperatures of ions and electrons. We also 
add neutral and ionized molecules. Molecules are more 
abundant in the Enceladus torus where low plasma density 
allows H2O to escape from Enceladus largely intact, whereas 
the more energetic local plasma interaction at Io results in 
dissociation of SO2 [Dols et al, 2008]. Moreover, thermal 
electrons 2 eV) throughout the Enceladus torus dissoci- 
ate H2O approximately ten times less easily than thermal 
electrons (« 5 eV) in the Io torus dissociate SO2 ( V. Dols, 
personal comm.). 

Previous models of molecular chemistry in Saturn's mag- 
netosphere were driven by Pioneer 11 and Voyager 1 and 2 
observations [Frank et al., 1980; Trainor et al., 1980; Wolfe 
et al, 1980; Bridge et al, 1981, 1982; Sittler et al, 1983]. 
Richardson et al. [1986] showed the importance of recom- 
bination under conditions of slow radial transport. Using 
essentially the same chemical reactions as Richardson et al. 
[1986], Richardson et al. [1998] determined ion and neu- 
tral lifetimes within Saturn's inner magnetosphere (< 12i?s; 
Rs = Saturn radius = 6.0 x 10 9 cm), constrained by HST ob- 
servations of the extended OH cloud (see their Figure 2, and 
references therein). They solved the rate equations for num- 
ber densities while also solving the radial diffusion equation, 
but energy conservation was not considered. Jurac et al. 
[2002] and Jurac and Richardson [2005] further improved 
on these models by considering neutral cloud expansion, and 
solved for plasma and neutral distributions self-consistently 
in order to study the source of water within Saturn's inner 
magnetosphere . 

Our model concentrates on the molecular chemistry. We 
start with a uniform box and characterize transport by just 
a time scale. However, we do consider energy balance. More 
importantly, unlike the above models, we retain HaO + in our 
model, which proves to be a significant component. 

The purpose of this paper is to summarize the sensitiv- 
ity of the chemistry of Enceladus's torus to several parame- 
ters. The parameters we investigated are hot-electron tem- 
perature, hot-electron density, H2O source rate, ion radial- 
diffusion time scale, and proton temperature anisotropy. Be- 
ing lighter, protons are less bound to the equator [Bage- 
nal et al, 1980; Sittler et al, 2008; Persoon et al, 2009]. 
This means protons spend only a fraction of the time in- 
teracting with the heavy ions and molecules. This effect is 
simulated with a 'proton dilution' factor (Section 3). We 
search for values of these parameters leading to thermal- 
electron temperature, thermal-electron density, and water- 
group (W+ ee 0+ + OH+ + H 2 0+ + H 3 0+) ion-to-proton 
ratio consistent with available Cassini data (Section 2.2). 

The observations used to constrain our model and to de- 
fine the parameter space are given in Section 2. The model 
is described in Section 3. The best fit (baseline solution) 
and the procedure used to find it are discussed in Section 
4. Model sensitivity to each of the parameters is discussed 
in Section 5. Finally, the importance of hot electrons with 
regard to water-group ion composition is demonstrated in 
Section 6. 

2. Observations 
2.1. Parameters 

Here we present the observations used to bound the pa- 
rameter search (Section 3). The baseline parameter values 
(listed in Table 1) are mentioned throughout this section 
and are discussed at length in Section 4. 

Neutral Source (iV src ): We have set the baseline 
source of water to be iV S cr = 2.0 x 10~ 4 molecules cm~ 3 s _1 . 



A torus centered on Enceladus with cross section (2Rs) 2 
has a volume of ^ 2n(4R s )(2R s ) 2 = 2.2 x 10 31 cm 3 , giv- 
ing a volumetric source rate of 4.4 x 10 27 H2O molecules 
s , or 130 kg s _1 . If one chooses a smaller or larger 
torus volume, the net neutral source rate is adjusted ac- 
cordingly. The best estimates of the Enceladus source come 
from Cassini observations. Hansen et al. [2006] estimated 
150 < iV src /(kg s _1 ) < 350 from two stellar occultation ob- 
servations of the Enceladus plume with Cassini UVIS. Tokar 
et al. [2006] inferred a source rate of 100 kg s _1 from the 
CAssini Plasma Spectrometer (CAPS) data acquired during 
the E2 Cassini flyby of Enceladus on 14 July 2005. 

Earlier source estimates come from models of the neutral 
clouds. Jurac et al. [2002] included the effects of collisional 
heating and developed a model to simulate the morphology 
of the extended OH cloud. They used neutral lifetimes de- 
rived from a two-dimensional model by Richardson et al. 
[1998] to determine the neutral H2O source responsible for 
the OH radial profile constrained by 1996 HST faint-object 
spectrograph observations. They found a total water source 
required to maintain the OH cloud of N src — 112 kg s _1 , 93 
kg s _1 coming from the orbit of Enceladus. 

In a later paper, Jurac and Richardson [2005] improved 
their model, treating plasma and neutrals self-consistently 
by tracking neutrals with a Monte Carlo algorithm and 
transporting plasma diffusively. They found a total water 
source rate of « 300 kg s _1 . A similar result was found by 
Burger et al [2007] with a three-dimensional Monte Carlo 
simulation of neutrals constructed to simultaneously model 
the Ion and Neutral Mass Spectrometer (INMS) and UVIS 
observations made during the E2 Cassini flyby of Enceladus. 

Radial Transport Time Scale (T"t rans ): We have 
set the baseline radial transport time scale to be 
n b rtn S linc = 26 days. Richardson et al. [1998] estimated 
Ttrans ~ 23 days from Voyager-era data and HST OH ob- 
servations. Using the radial velocities of Sittler et al. [2006] 
(their Figure 9), one finds a diffusion time scale at Ence- 
ladus of w 2Rs/vn(r = 4i?s) ~ 12 days. This more rapid 
diffusion might suggest that the magnetosphere was com- 
pressed on Saturn Orbital Insertion (SOI), which would in- 
crease the angular speed to beyond corotation as momentum 
is conserved, ultimately resulting in enhanced radial outflow 
velocities (reduced transport times). Thus, the Sittler et al. 
[2006] radial velocities based on SOI data may not represent 
the entire magnetosphere, and not for all epochs. Radial 
convection may also be superimposed on the diffusive mo- 
tions throughout Saturn's inner magnetosphere due to flux 
tube interchange instabilities [Rymer et al, 2008]. A major 
goal of this study is to explore the consequences of such a 
wide range of time scales for radial transport. 

Hot-electron Temperature and Fraction (T e h, /eh): 
Our model gives a baseline value for the temperature of 
hot electrons of T c b h asclinc = 160 eV. The Cassini Radio and 
Plasma Wave Science (RPWS) measurements by Moncu- 
quet et al. [2005] indicated a hot-electron component within 
3-57? s with 40 < T ch /(eV) < 90. From Cassini CAPS 
ELectron Spectrometer (ELS) observations acquired dur- 
ing the SOI period, Young et al. [2005] found that T oh 
ranged from « 500-1000 eV inbound and w 800-3000 eV 
outbound (SOI), indicating a strong longitudinal and/or 
temporal dependence. Lewis et al. [2008] performed mo- 
ment calculations on CAPS-ELS data with two differ- 
ent methods and found 300 < T ch /(eV) < 2000 eV and 
1500 < T ch /(eV) < 4000. Schippers et al. [2008] combined 
CAPS-ELS and Cassini Magnetospheric IMaging Instrue- 
ment (MIMI) data and found significant scatter in T eh and 
calculated 200 < T h/(eV) < 2000 at the closest approach 
of5.4fl s . 
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Our baseline value for the fraction of the electron den- 
sity in the hot component is found to be f^ schnc = 0.46%. 
Young et al. [2005] found that / e h ranged between 0.01% 
and 5% within 3-5-Rs- As with T e h, / c h varied significantly 
between the inbound and the outbound data. Lewis et al. 
[2008] found / e h from their 3-D moment calculation to be 
< 1%. Schippers et al. [2008] calculated f ch as low as 0.1% 
(inbound) and as high as 0.3% (outbound). 

At the time of Sittler et al. [2008], the authors felt 
that hot-electron parameters measured by CAPS ELS were 
highly uncertain in the vicinity of Enceladus, due to pene- 
trating radiation. Instead, they used the Moncuquet et al. 
[2005] RPWS observations, which found that T ch w 50 eV. 
For the hot-electron density, they used n c h ~ 0.1 cm~ 3 from 
the Sittler et al. [1983] Voyager observations (which are not 
affected by penetrating radiation). Sittler et al. [2008] com- 
bined the Moncuquet et al. [2005] Cassini RPWS data for 
T e h and the Voyager data for n c h to compute the total (effec- 
tive) electron temperatures T c . Since Voyager and Cassini 
SOI were so different in time, they used only the total elec- 
tron temperatures for their reaction rates. 

Proton Dilution (/h+ ) '■ Previous work by Sittler et al. 
[2008] (their Figure 4) showed that « 2/3 of the proton pop- 
ulation are distributed within a distance of IRs from the 
centrifugal equator at the orbit of Enceladus. Protons are 
pulled above the equator by the ambipolar electric field and 
do not couple efficiently to the heavy water-group ions. 

In Section 5.1.1, we show that the proton abundance is 
strongly coupled to the hot-electron population via impact 
ionization. Our model is consistent with any value of the 
proton dilution parameter / H + between 0.7 and 1.0. Values 
of / H + = 1 or represents the cases where no or all newly 
created protons are excluded from the model. Thus, to sim- 
plify the present analysis, jbaseime has been set to unity. 

2.2. Constraints (n e , T e , W+/H+) 

Parameter combinations are evaluated by comparing the 
corresponding model output to the following observations. 
We choose n c = 60 cm" 3 , T c = 2eV, and W+/H+ = 12 as 
the initial model constraints (Section 4). 

Total Electron Density (n e ): Gurnett et al. [2005] re- 
ported 20 < n /(cm~ 3 ) < 100 within 3-57?s during the 
approach and first orbit around Saturn (SOI) from the 
Upper Hybrid resonance Frequency (UHF), acquired by 
the RPWS instrument. Moncuquet et al. [2005] consid- 
ered RPWS Quasi- Thermal Noise (QTN) on SOI and found 
40 < n c /(cm- 3 ) < 70 within 3-5/fe. Persoon et al. [2005] 
used the RPWS UHF to determine n c from five later orbits. 
They discovered variability inside ~ 5Rs, with n c ranging 
from 35-105 cm -3 . In a later paper, Persoon et al. [2009] 
developed a diffusive equilibrium model from RPWS and 
CAPS data acquired on 50 passes through Saturn's inner 
magnetosphere from 30 June 2004 to 30 September 2007, 
and calculated an electron density of over 50 cm -3 near the 
orbit of Enceladus. 

Lewis et al. [2008] derived 6 < n /(cm~ 3 ) < 20 within 
3-5i?s from the moment calculations with CAPS-ELS elec- 
tron distribution data. Their analysis was based on SOI 
data when CAPS was not fully actuating [Sittler et al., 
2006, 2007]. Also, Cassini is a three-axis stabilized space- 
craft with a fixed field of view. Hence, sampling is lim- 
ited, and moment calculations are not as straight forward 
as they are for spinning spacecraft, as discussed by Lewis 
et al. [2008]. Additionally, Cassini was likely negatively 
charged on SOI [Young et al, 2005], resulting in a lower- 
than-expected n e from the moment calculations. 

Schippers et al. [2008] performed a multi-instrument anal- 
ysis of the electron populations for several orbits and found 
n c Ri 10 cm~ 3 from their own CAPS-ELS analysis, but ad- 
mitted that a negative spacecraft potential inside 9Rs likely 
resulted in an underestimate of the thermal-electron density. 



Instead they used the RPWS UHF analysis by Gurnett et al. 
[2004] in this region, where n c « 50 cm -3 near the orbit of 
Enceladus. 

From consideration of the above-cited observations, we 
take a value of n c — 60 cm~ 3 for the total electron density 
in our model. 

Thermal-electron Temperature (T e ): Moncuquet 
et al. [2005] derived 1 < T e /(eV) < 4 within 3-57? s horn 
RPWS data. The SOI CAPS analysis by Young et al. [2005] 
found that 2 < T e /(eV) < 20. By comparison, the CAPS 
analysis by Lewis et al. [2008] found that T e ranged from 
roughly 2 to 4 eV within 3-57?s using one method, and be- 
tween 1 and 4 eV using another. Sittler et al. [2006, 2007] 
found T e ~ 1.5 eV between 3.5 and 4.5 Rs. Schippers et al. 
[2008] calculated T c « 2 eV at closest approach (5.4i?s) 
from a multi-instrument analysis. This value of T e = 2 eV 
is the one we use for the thermal-electron temperature in 
our model. 

Ion Composition (W+/H+): Sittler et al. [2008] cal- 
culated H + and W + densities from the CAPS SOI results 
of Sittler et al. [2006, 2007]. Their Figures 5 and 6 show 
that W + /H + ranged from £s 3-15 over the torus (3-5i?s)- 
The CAPS SOI analysis of Young et al. [2005] found 6 < 
W + /H + < 35 within 3-5ife- The Young et al. [2005] anal- 
ysis was based on non-coincident CAPS IMS singles data, 
whereas the Sittler et al. [2008] results were based on co- 
incident time-of-flight data. It is challenging to distinguish 
protons from water-group ions in the singles data, while pro- 
tons are well-separated from the water-group ions in their 
respective time-of-flight channels. Observations based on 
IMS singles data are thus likely to overestimate the W + /H + 
ratio. 

Wilson et al. [2008] presented a forward modeling tech- 
nique to CAPS data for dayside equatorial orbits between 
5.5 and 117?s to calculate ion densities and temperatures. 
Extrapolating their results down to 5Rs gives W + /H + ~ 15. 
Persoon et al. [2009] used anisotropy measurements from 
CAPS and temperature measurements from RPWS and 
CAPS to estimate equilibrium distributions within Saturn's 
inner magnetosphere. Near the Enceladus torus, they found 
W+/H+ « 10. We choose a value of W+/H+ = 12 for the 
water-group to proton density ratio in our model. 

3. Model 

Here we present a model based on the Neutral Cloud 
Theory (NCT) model described in Delamere and Bagenal 
[2003], which was developed to address the variability of 
plasma conditions in the Io torus. The Delamere model 
was based on the earlier NCT models of Shemansky [1988] , 
Barbosa [1994], Schreier et al. [1998], and Lichtenberg et al. 
[2001]. The tools developed by Delamere and Bagenal [2003] 
to study sensitivity of the plasma-dominated environment at 
Io are utilized here to study the neutral-dominated water- 
based Enceladus torus. 

The model is 0-dimensional and homogeneous. In this 
paper, we calculate steady-state densities and temperatures 
of ions and neutrals originating from a pure H2O source. 
When volumetric quantities are reported, we have adopted 
the volume used by Delamere et al. [2007] [2tt(47?s)(2-Rs) 2 « 
2 x 10 31 cm 3 ], which roughly corresponds to a torus of minor 
radius life centered on Enceladus's orbit at 4Rs- Because 
pickup temperatures vary with radial distance along Sat- 
urn's equatorial plane (appendix, Eqs. 6 and 7), the scaling 
is only approximately valid for the span considered here of 
3-5i? s . 

The basic equations [Barbosa et al, 1983] for number den- 
sity and energy density for species a are 
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and 



d(ln a T a ) 
at 

The S a 's and £ Q 's represent source rates and loss rates, 
respectively, for species a. Following the convention of De- 
lamere and Bagenal [2003], the factor of 3/2 in Eq. 2 will be 
dropped henceforth, so that we are actually solving for an 
'effective' temperature rather than energy. The temperature 
is described as effective because for pickup ions, the temper- 
atures perpendicular to the magnetic field are expected to be 
greater than the parallel temperature (for observations and 
further discussion see Richardson and Sittler [1990]; Mon- 
cuquet et al. [2005]; Sittler et al. [2006, 2007]; Tokar et al. 
[2008]; Per soon et al. [2009]). A complete discussion of Eqs. 
1 and 2 for ions, electrons, and neutrals can be found in the 
appendix. 

The rate at which the H 2 particles are introduced into 
the model (N BTC ) is a free parameter. Chemical pathways en- 
compass a set of reactions involving H, H + , H2, , O, + , 
0++, OH, OH+, H 2 0, H 2 + , and H 3 0+ (see the appendix 
for a complete list of reactions). Note that we assume 50% 
of the hydrogen produced from impact dissociation of H2O 
and OH has enough energy to escape our model [Richardson 
et al., 1998]. Neutrals are assumed to be cold, having only 
bulk motion. They are not collisionally-heated in the model, 
and here we assume that neutrals created from ion charge 
exchange have velocities greater than the escape speed from 
Saturn and are ejected from the model. Steady-state num- 
ber and energy densities are found for each species by solving 
Eqs. 1 and 2 iteratively using a modified Euler method with 
second-order accuracy. 

At the heart of this work is a sensitivity investigation of 
model output within the parameter space spanned by: 

— Neutral source rate (iV src /10~ 4 cm~ 3 s~^): 0.2 -^3.0 

— Hot-electron temperature (T eh /eV): 20 — » 400 

— Hot-electron fraction (/ c h = n c h/n c ): 0.05 — > 1.0% 

— Radial transport time scale (r tr ans/days): 2 -^60 

— Proton dilution factor (/h+) : 0.7 — > 1.0. 
These ranges reflect the broad set of observations given in 
Section 2.1. The proton dilution factor, / H +, removes pro- 
tons from the model, and has been implemented by modify- 
ing the source term in Eq. 1 for protons: 



dn 



11+ 



dt 



— /h+<Siti,H+ ~£ m ,H+- 



In reality, the heavy ion abundance peaks near the centrifu- 
gal equator, while the proton abundance peaks well away 
from the equator and out of our model domain [Persoon 
et al, 2009]. We apply the above equation to crudely ad- 
dress and investigate this phenomenon. 

4. Baseline Solution 
4.1. Procedure 

Initially we set / H + = 1 (see further discussion in Sec- 
tion 5.1.1) and explore the space spanned by / c h, T e h, Ttrans, 
and JVsrc alone. The model was run with a random-walk 
Monte Carlo algorithm to find the parameter space coordi- 
nates yielding the best agreement between model output 
and the constraints. The following have been chosen as 
the constraints on the model (Section 2.2): n e = 60cm~ 3 , 
T c = 2eV, W+/H+ = 12. In Section 5.1.2 we accommodate 
a wider range of observations and investigate how composi- 
tion is affected by these choices of n c , T e , W + /H + . 

We define best agreement as the smallest total fractional 
difference between the model output and the constraints: 



Model; 



Constraint; 



where i — n c , T c , W + /H + in the present case. The baseline 
solution (parameter combination with the smallest /diff) was 

(2) found by starting the random-walk algorithm from a point 
in parameter space near the global minimum in /diff- The 
procedure used to find the global minimum is discussed in 
Section 5. Model output was evaluated and a step was ran- 
domly taken in the direction of one of the four parameters. 
The step sizes were also random in length and constraint- 
dependent. For example, the step size for the hot-electron 
fraction ranged from 0-0.001, while the step size for the 
transport time ranged from 0-86,400 seconds (1 day). Every 
value throughout the step intervals had equal weight. The 
model output was then evaluated (/diff calculated), and the 
procedure was repeated until a minimum in /diff was found. 
The random walk led to the baseline parameter combina- 
tion, where faT 1 ™ w 0. 

The solution and corresponding model output are given 
in Table 1. Lists of the values of lifetimes of each species 
controlled by the primary source/losses mechanisms are pre- 
sented in Table 2, and lifetimes listed by each separate re- 
action can be found in the appendix. 

4.2. Results 

Table 1 presents the model output for the densities and 
temperatures of all species. We find a torus composition that 
is dominated by neutral species with roughly equal amounts 
of H, O and OH (~ 700 cm~ 3 each) with a lesser amount of 
water molecules (190 cm ) and trace amounts of H2. The 
ion species are dominated by H 2 0+ and OH+ (^20 cm" 3 ) 
followed by H 3 0+ and 0+ (w 10 cm" 3 ), H+ (w 5 cm" 3 ) and 
trace amounts of ++ and Hj. 

When we compare our baseline OH density (770 cm -3 ) to 
Figure 3 of Jurac et al. [2002] wc find very similar values. 
Their model was developed to simulate the morphology of 
Saturn's extended OH cloud, as measured by HST, October 
2002. They found an OH density of > 750cm~ 3 centered 
on the orbit of Enceladus (see also Jurac and Richardson 
[2005]). This is an independent test of our results since no 
radiative constraints were used to determine the baseline 
solution. 

(3) Our model shows that all ion species have temperatures 
close to their initial pickup temperature, consistent with 
negligible loss of energy via radiation and Coulomb colli- 
sions with electrons. The thermal coupling time between 
electrons and ions derived from the model is « 60 days. Be- 
cause ions are transported out of the box in 26 days, they do 
not efficiently transfer energy to the electrons. Oxygen ions 
are picked up by the corotating magnetosphere with a tem- 
perature of 38.4eV (appendix, Eq. 7). OH, H 2 0, and H3O 
ions are picked up with 40.8, 43.2, and 45.6 eV, respectively. 
As shown in Table 1, very little of these heavy ions' thermal 
energy has been transferred to the thermal electrons. This 
result has also been established in Figure 2 by the small 
energy coupling (2.4%) between ions and electrons. 

The water-group temperatures from our model (« 38- 
42 eV) are warmer than the CAPS data suggest. Sittler 
et al. [2006, 2007] observed that T ± w+ « 35-40 eV near the 
Encaladus torus. With their anisotropy of (Tx/Ty) w + ~5, 
the effective water-group temperature is reduced to T w + = 
(2T± + Tj| )/3 fa 27 eV. The discrepancy between our model 
and the data may be explained by a sub-corotating plasma 
torus near Enceladus' L shell. An w 20% sub-corotation of 
the plasma flow, as measured by CAPS at 4Rs (R. Wilson, 

(4) personal comm.), would reduce pickup energies and may ac- 
count for the difference between our model temperatures 
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and the Sittler et al. [2006, 2007] observations (40 eV and 
27 eV, respectively). 

In our model, we have assumed that the ion velocity dis- 
tributions are isotropic and thus cannot comment on the 
parallel and perpendicular temperatures individually (Sec- 
tion 3). Because the data suggest that the water-group 
ion has an anisotropy of T±/Tu « 5 and the protons have 
an anisotropy of T±/Tu ~ 2 [Richardson and Sittler, 1990; 
Moncuquet et al., 2005; Sittler et al., 2006, 2007], we hope to 
include anisotropic ion velocity distributions in the future. 

The flow of mass and energy is shown in Figures 1 and 2 
(contributions from individual species are listed in the ap- 
pendix). Mass is introduced into the model by way of H2O 
only and leaves when ions radially diffuse, charge exchange 
with neutrals, or recombine with electrons (for singly-ionized 
species). Recombination and charge exchange represent 
mass sinks because the ions become 'fast neutrals,' assumed 
to possess enough velocity to escape the torus. We find that 
94% of the particles leave the torus as fast neutrals and 6% 
as diffusing ions. 

Energy is introduced almost entirely by pickup ions, 
and a small amount (1.8%) comes from Coulomb collisions 
between the hot electrons and the thermal electrons/ion 
species. Pickup ions represent an energy source due to the 
velocity difference between neutrals and ions in the Ence- 
ladus torus; a freshly-ionized neutral is accelerated to mag- 
netospheric corotation via current systems established be- 
tween Saturn's ionosphere and the Enceladus torus. Fresh 
pickup ions may be produced either by electron impact or 
charge exchange between an ion and a neutral. Energy is 
carried away in the model by the fast neutrals, diffused ions, 
and radiation induced by electron-impact excitation. Figure 
2 indicates that most energy (83%) leaves the model with 
the fast neutrals. 

Figures 1 and 2 can be compared to the energy and par- 
ticle flow diagrams produced by the oxygen-based model of 
Delamere et al. [2007] (their Figure 1). In the partitioning 
of particle flow between radial transport and fast neutrals, 
their simplified model is remarkably similar to our solution, 
though their neutral source rate (4 x 10~ 4 cm~ 3 s _1 ) is twice 
as strong. 

The total energy flowing through our model is roughly 
40% of that found by Delamere et al. [2007] (9.4eVcm -3 s" 1 
compared to 23eVcm~ 3 s _1 ). They used a smaller hot- 
electron fraction (/ e h = 0.3%) and a higher hot-electron 
temperature (T e h = 1000 eV) than the values we used to 
produce Figures 1 and 2. In addition, their transport time 
was considerably longer at 45 days. 

Regarding energy output, we find that more energy is 
transported out of the torus by ions than by fast neutrals 
when compared to Delamere et al. [2007], though we agree 
on the radiated power partition of a few percent. Energy in- 
put is remarkably different because we have included molec- 
ular chemistry, while Delamere et al. [2007] included atomic 
oxygen only. Hot-electron thermal coupling plays a bigger 
role as an energy source with 1.8% of the total energy input 
compared to 0.5% in Delamere et al. [2007]. In our model, 
charge exchange only marginally exceeds the combination of 
photoionization and impact ionization as a means of adding 
fresh pickup ions to the system, whereas in Delamere et al. 
[2007] charge exchange was found to dominate these mecha- 
nisms by a factor of 19. One reason for this major difference 
is the relative ease at which OH is ionized via electron impact 
with respect to O. This reaction can contribute greatly to 
the overall energy budget since OH is the dominant species 
(Table 1). Photoionization of OH and O occurs at roughly 
the same, much lower rate for the baseline case. If oxy- 
gen were the only species in the model, energy input would 
be a competition between the highly-likely resonant charge 
exchange between O and + and the order-of-magnitude- 
less-likely photo- plus impact ionization of O. The reaction 
rates supporting the above argument can be found in the 
appendix. 



Dissociative recombination represents an important 
plasma sink to the Enceladus torus. This process has a pro- 
found effect on the neutral-to-ion ratio in the torus. The De- 
lamere et al. [2007] model found 12, but they 
argued that this ratio represents a lower limit since their 
model is oxygen-based and does not include dissociative- 
recombination reactions; the recombining ions' neutral prod- 
ucts have escape velocities and leave the model, just as re- 
combining atomic ion species do. We find n ncu t/ni ons = 40 
with the full water-based molecular chemistry. 

Sittler et al. [2008] showed that the combination of cold 
electrons (^ 1 eV) and the dominance of molecular ions 
over atomic oxygen ions near Enceladus [Young et al, 2005] 
drives a rapid dissociative-recombination time scale. They 
also showed that this process increases the neutral-to-ion 
production ratio to « 50 near Enceladus (their Figures 18 
and 20). 

To illustrate the importance of dissociative recombination 
in our model, we increased T e by increasing the pickup tem- 
perature in the Enceladus torus (appendix, Eq. 7). Indeed, 
if the parameters are fixed at the baseline values (Table 1), 
and T c is increased to 6 eV (as in the ion-dominated Io 
torus), dissociative recombination continues to prevent ions 
from dominating neutrals (Section 4.2.2). At Saturn, im- 
pact ionization by thermal electrons cannot compete with 
dissociative recombination as an ion sink, even when T c = 6 
eV. 

Derived Quantities: The field strength in the Ence- 
ladus torus is 325 nT [Dougherty et al., 2006], giving a 
plasma beta of 

B 2 - /on 

where the summation is taken over all charged species, in- 
cluding both the thermal- and hot-electron populations (Ta- 
ble 1). This plasma beta is consistent with the Sittler et al. 
[2008] analysis, which found f3 between 3 and 5Rs to be 
between 0.1 and 5%. The Alfven speed is given by 

v A = B/x/4^o = 230kms~\ (6) 

where p — X^j=ions n o m o- If the plasma is at full corotation 
at the orbit of Enceladus [as has been assumed in the model 
for calculating E pu (appendix, Eq. 7)], then the Alfven Mach 
number is 

M A = v^/va = (4Rs)n s /v A = 0.17. (7) 

This should be compared to Sittler et al. [2008], who find 
0.01 < M A < 0.5 between 3 and 5R S . Sittler et al. [2008] 
use ion-electron fluid parameters as boundary conditions 
to solve for ion field-line distributions throughout Saturn's 
inner-magnetosphere. The fluid parameters are derived from 
CAPS data acquired during the approach phase of the SOI 
period [Sittler et al, 2006, 2007]. In a future study, we 
will solve the radial-transport equation and present a self- 
consistent map of the ion distribution throughout Saturn's 
inner magnetosphere (as performed by Richardson et al. 
[1998]). 

4.2.1. Lifetimes 

The lifetimes are listed by mechanism in Table 2. The 
lifetimes are listed by reaction in the appendix. The diffu- 
sion lifetime (r tralls ) has not been well-constrained by this 
study. In Section 5, we find that the model is consistent 
with a diffusion time scale of 12 days and longer. Because 
our study cannot place an upper limit on r tra ns, it is not 
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possible to say which mechanisms occur more rapidly than 
radial transport and are therefore more important. 

Our model does not account for collisional heating mech- 
anisms that would give neutrals enough speed to escape the 
model domain. Thus, we have assumed that the time scales 
for such neutral escape are longer than the time scales for 
the included mechanisms. Based on the Enceladus torus 
study by Farmer [2009], wc now examine this assumption 
and conclude that the effects of collisional heating should be 
included in future iterations of our model, especially in the 
case of OH. 

Neutral H2O has an effective cross-section to neutral- 
neutral collisions due to dipole-dipole interactions such that 
an H2O molecule will be transported outside our modeled 
torus (>57?s) after & 40 days (Figure 2 in Farmer [2009]). 
Then, according to Table 2, only impact dissociation and 
charge exchange are important loss mechanisms for H2O. 
A similar time scale (40 days) may also limit OH lifetimes 
because its induced dipole is comparable to that of H2O. Un- 
der this assumption, collisional heating would be the most 
important loss mechanism for OH. 

The Farmer [2009] result can also be used to estimate 
oxygen lifetimes against neutral-neutral collisions. Be- 
cause the geometric cross section for H2O is a factor of 
10 smaller than the induced-dipole cross section (o"h 2 o = 
5 A 2 -><TH 2 d = 54 A 2 , Farmer [2009]), and because the col- 
lision frequency is proportional to a, one might expect that 
O will take roughly 10 times longer to be scattered outside 
the torus. That is, only those mechanisms occurring on a 
time scale of < 400 days (charge exchange and impact ion- 
ization) would occur before neutral collisions remove oxygen 
from the torus. 

4.2.2. Dominant Chemistry 

Reactions occurring more frequently than 10~ 6 cm~ 3 s _1 
, shown in Table 3, are the ones of primary importance to the 
torus chemistry. The appendix lists the full set of reactions 
and reaction rates for the baseline case. For example, if we 
run the model with every reaction in the appendix turned 
on (and the parameters set at the baseline values) all densi- 
ties and ion temperatures are within 3% of the results from 
the calculation using the refined set presented in Table 3. 

The most dominant reaction is impact dissociation of 
H2O by hot electrons (H2O + e h — > OH + H + e). Dissocia- 
tion of H2O by hot electrons occurs so frequently because of 
the relatively large reaction rate at T c b h asolmo = 160 eV (ap- 
pendix) as well as the high baseline density of H2O (Table 
1). Several other charge-exchange, photolytic, and electron- 
impact reactions are competitive behind H2O impact disso- 
ciation. 

Impact ionization by hot electrons contributes roughly 
the same amount of energy (22% of total input) via mag- 
netospheric pickup as does photoionization (16%). Im- 
pact ionization by thermal electrons is a minor source of 
energy to the torus (~ 1%). In fact, only one reac- 
tion in Table 3 involves thermal-electron impact ionization 
(OH + e^ OH+ + 2e). Delamere et al. [2007] find that 
charge exchange is far more important than photo- and im- 
pact ionization combined as a torus energy source. Here we 
find that charge exchange is only marginally more impor- 
tant (ss 60%) than the combination of these other ionization 
sources (ss 40%) for providing fresh pickup ions to the torus. 
This discrepancy (discussed in Section 4.2) is largely due to 
the fact that the earlier model did not include molecular 
chemistry. 

4.3. Sensitivity 

Contour plots of the total fractional difference (/diff, Eq. 
4) between the model output and the constraints have been 
created for every parameter combination and are shown in 
Figure 3. The intersection of the dashed lines indicates the 



baseline solution (Table 1). In each case, the remaining 
three parameters are fixed at the baseline values. The shad- 
ing inside the /diff = 1 contour is intended to guide the eye 
when comparing one panel to another. Because the remain- 
ing two parameters are fixed in each panel, these plots show 
the sensitivity of /diff to each parameter individually. 

The source rate is inversely related to rtrans (Panel 1) 
and /eh (Panel 2). The trend between iV src and Ttrans Can 
be understood as a balance between source and sink; plasma 
taking longer to diffuse out of the model must be accompa- 
nied by a decrease in H2O. Similarly, an increase in hot elec- 
trons results in higher ionization. A smaller source rate is 
required to maintain the torus composition consistent with 
the constraints on n e and W + /H + . 

The total fractional difference, /diff, strongly depends on 
/eh, with a strongly pronounced valley in all cases (Panels 
2,4,5). The hot-electron population is critical for ionizing 
H efficiently to obtain W + /H + w 12 (and hence, minimizing 
/diff). The hot electrons are equally necessary for attaining 
a higher overall ionized composition, thereby increasing n e , 
and for heating the thermal electrons (T e ) via Coulomb cou- 
pling. 

A similar dependence exists for T c h, except that the 
strong dependence is at low T c h only (Panels 3,5,6). Be- 
yond T c b h asolino (160 eV), /diff is roughly independent of T ch . 
The other three parameters dominate variation in f di g when 
the hot-electron population is sufficiently energetic (> 160 
eV). 

Delamere et al. [2007] present sensitivity contour plots 
of the neutral-to-ion ratio and thermal-electron tempera- 
ture from their oxygen-based model (their Figure 2). We 
have generated similar contours (not shown) and find that 
the torus composition is neutral-dominated despite T e ap- 
proaching 6 eV, whereas Delamere et al. [2007] find that 
ions dominate when T c is as low as 3 eV. The key differ- 
ence between their model and ours is that we have included 
a complete water-based, molecular set of reactions, and we 
have included the effects of dissociative recombination not 
present in the Delamere model. 

Because we are varying only two parameters at a time in 
each panel of Figure 3, we are not necessarily finding the 
highest quality of fit (smallest /diff) throughout. This con- 
cern is addressed in Section 5. 

5. Grid Search 

In search of the best possible match to our three observ- 
able constraints, we have explored the full 5-dimensional 
parameter space (/ eh , T c h, / H +, N SIC , r tra n S ) by solving Eqs. 
1 and 2 systematically for a large set of parameter combina- 
tions. We divided the domain for each parameter (Section 
3) into 30 discrete, uniformly-spaced values. We created a 
table of model output values of densities and temperatures 
corresponding to every parameter combination. This table, 
or grid, is an extremely useful tool for comparing output 
from our model to any combination of observations (Section 
5.1.2). The computationally-expensive procedure of creat- 
ing the grid will only need to be repeated when additional 
or updated chemical reactions are introduced to the model. 
The downside of the grid search is limited resolution, and 
to double the resolution would require 2 5 times more com- 
putational time. Fortunately, the calculation can be done 
in parallel, and wall-clock time can be reduced linearly with 
the number of computer processors. 

5.1. Results 

In the present case, we have used the grid to generate sen- 
sitivity contours for every parameter combination by search- 
ing for the smallest /diff everywhere using the constraints in 
Section 4: n c = 60 cm~ 3 , T c = 2eV, W+/H+ = 12. The pri- 
mary purpose of these contours is to establish the following 
parameter limits consistent with the primary constraints: 
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neutral oxygen emission lines. If measured, Pijv would be a 
powerful constraint for our model [Esposito et al, 2005]. 

Radial transport, or mass loading (M, Panel 13), is de- 
fined here as the total radially-transported mass: 



1-5 < AW(10- 4 cm" 3 s- 1 ) < 2.7 
12 days < r tra n S 
0-3 < /eh(%) < 0.9 
100 < T ch /eV < 250 
/h+ < 1- 

We will also demonstrate that the solution space found by 
searching the grid is smoothly varying and well-behaved. 

In Figure 4, we present the sensitivity contours between 
the parameters /ch and T en . Panel 1 is the contour plot for 
the smallest obtainable total fractional difference within the 
domain defined in Section 3. Recall that /diff is the sum 
of fractional differences between model output and primary 
model constraints for n c , T , and W + /H + (Section 2). To 
illustrate the individual contributions to /diff, we show T e , 
n c , and W + /H + in Panels 2, 3, and 4 respectively; the dot- 
ted contours define the primary constraints used to evaluate 
/diff. The values of the parameters r tr ans, N SIC , and / H + 
corresponding to the model output for n c , T , and W + /H + 
are plotted in Panels 5-7. We also present the contours for 
model output quantities: water-group composition (Panels 
8-11), total UV power (Panel 12), and mass radial transport 
rate (Panel 13). 

The /diff = 0.05 contour in Panel 1 has been shaded in 
gray and over-plotted in Panels 5-13 to indicate a parameter 
subspace consistent with the observations. This contour can 
be expanded or reduced to reflect observational uncertainty. 
The random-walk algorithm used in Section 4 to find the 
baseline solution was initialized with parameter coordinates 
near the global minimum inside this solution space. The 
baseline solution has been indicated throughout Figure 4 by 
the intersection of dashed lines. Though the baseline solu- 
tion is encompassed in the solution space, it is not unique. 
The solution space is bounded by a range of the parameters 
for which a solution can be found that agrees with the data 
approximately as well as the baseline solution discussed in 
Section 4. From Figure 4, we find that 0.3 < / oh (%) < 0.9 
and 100 < T oh /(cV) < 250. 

Because parameter combinations are evaluated according 
to model output for n e , T c , and W + /H + alone, one can- 
not rule out a priori that adjacent points in these contours 
sample wildly different values of r trans , iV src , and / H +. That 
the contours for these constraints (Panels 5-7) and for the 
composition (Panels 8-13) are smoothly varying and well- 
behaved proves that this is not the case. 

Because the different water-group ion species have simi- 
lar masses and because composition varies significantly with 
longitude in the Enceladus region [Williams et al., 2008], 
the water-group density ratios (Panels 8-11) are difficult 
to distinguish in the CAPS observations. However, Sittler 
et al. [2008] have shown, based on CAPS SOI data, that 
H 3 0+/W + ~ 0.45, H 2 0+/W + « 0.15, and OH+/W+ w 
0+/W+ « 0.2 near the orbit of Enceladus. Preliminary 
CAPS results by Williams et al. [2008] also suggest that 
H3O 4 " is the dominant water-group species in the Enceladus 
torus. 

INMS data acquired downstream of Enceladus suggest 
that the local chemistry, dominated by charge exchange, 
may represent a significant source of HsO + not included 
in our model [Cravens et al., 2009]. We are currently 
preparing a manuscript on the chemical interaction between 
Enceladus's water-based plumes and the corotating plasma 
in which we identify efficient chemical pathways by which 
H30 + is created. We hope to incorporate this correction 
into our model and to properly compare our results to the 
CAPS data in a future study. 

The UV power due to impact excitation (Panel 12) was 
calculated from line emissions and scales to roughly 1 GW 
for the entire torus. This estimate likely represents a lower 
limit since emission from H2O, for example, has not been 
considered in our model (see appendix). Most of the radi- 
ated power (« 80%) comes from the 1304, 1356, and 6300 A 



Radial Transport Rate = mjnj, (8) 

rtrans 

j— ion,e 

where a torus volume of 27r(4i?s)(2-Rs) 2 has been used, and 
rrij is the mass of ion species j. We find that ion transport 
may vary by a factor of six or more (« 8-50 kgs -1 ) and still 
be consistent with the solution space presented here. 

The solution space can be further constrained when bet- 
ter data for the parameters themselves [iV src , / c h, Teh, /h+ 1 r tr: 
(Section 2.1)] become available. Limits derived from such 
measurements could be used to define contour levels in Pan- 
els 5-7 that would limit the solution space in gray. The 
same can also be said for torus ultraviolet emission (Puv). 

The fractional difference contours in Panel 1 of Figure 4 
are much broader than in Panel 5 of Figure 3. The difference 
of course being that rtrans, N SIC , and / H + are fixed in Figure 
3 while in Figure 4 they are free. To help understand how 
the free parameters (Panels 5-7) have expanded the solution 
space, we offer the following interpretation. For clarity— and 
because the proton abundance is strongly coupled to hot- 
electrons (Section 5.1.1)— we restrict our discussion to trans- 
port time and source rate. Focus on the upper-right quad- 
rant of Panels 5, 6 and 7 in Figure 4. This is the region 
where the total fractional difference contours (Panel 1) have 
been broadened most markedly when compared to Figure 3. 
The transport time increases slightly (r trans f ) with increas- 
ing hot-electron temperature (T en T) an< i decreases sharply 
(rtrans I) as the hot-electron fraction increases (/ c h |). The 
neutral source rate is remarkably constant at a value near 
the baseline value of 2.0 x 10 -4 cm -3 s _1 , varying by only as 
much as 30% throughout the entire upper-right quadrant. 

The trends between /oh, and T eh are in part driven 

by the thermal-electron density. As the hot electrons in- 
crease in temperature beyond the baseline hot-electron tem- 
perature, 160 eV (T c h T)j the impact ionization rate of both 
hydrogen and water-group molecules decreases (Figure 8). 
To maintain the total ionization (n e ), the ions must remain 
in the torus longer (rtrans T)- Similarly, as the hot electrons 
increase in number (/ e h 1), both protons and water-group 
ions build up and must transport out of the model more 
rapidly (rtrans 1) to maintain n e . 

In Figure 5, we also include the same set of sensitivity 
contours, this time between r trans and / en . The solution 
space is bounded by 0.3 < / c h(%) < 0.9, but rtrans is only 
bounded from below at 12 days. There are three other sen- 
sitivity plots in the appendix with which to constrain rtrans, 
but rtrans is never constrained from above because recombi- 
nation and charge-exchange dominate the chemistry when 
rtrans ^ 26 days. That is, varying rtrans in this regime has 
no effect on the model output. The importance of recombi- 
nation has also been discussed in Sittler et al. [2008]. 

Sensitivity plots similar to Figures 4 and 5 for the eight 
remaining parameter combinations can be found in the ap- 
pendix. The /diff = 0.05 shaded contours in the appendix 
have been used to find the limits on the neutral source rate 
(iV src ) and the proton dilution factor (/ H +) presented above 
and in Section 7. 
5.1.1. Proton Dilution 

The grid search has revealed that / H + is a weak parameter 
in the model. That the model output is weakly dependent 
on / H + is illustrated in Figures 4 and 5. Panel 7 in each 
figure shows that small changes in Th and f ch can com- 
pensate for large changes in / H +. In fact, the full domain of 
/ H + (0.7-1.0) fits entirely into the shaded solution subspace. 
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Table 3 (and the reaction rates in the appendix) indicates 
that the dominant reaction for creating protons is impact 
ionization of hydrogen by hot electrons: 

H + c h -> H+ + 2e. (9) 

This reaction strongly couples / H + to / e h, reducing the pa- 
rameter space to four dimensions. The reaction H2O + eh — > 
H + + OH + 2e is only about 1/3 as effective as hot-electron 
impact ionization at producing fresh H + (appendix), while 
+ + H — > O + H + is independent of hot electrons alto- 
gether. 

Figure 6 illustrates that both the fit and the composi- 
tion depend weakly on the choice of / H + . The grid was 
searched to find the best combination of T e h, / e h, N STC , and 
Ttrans for three different values of / H +: 0.7, 0.85, 1.0. The 
combination yielding the best agreement is shown in the left 
panel— normalized to the baseline values in Table 1— and the 
model output for water-group composition is shown in the 
right panel. 

5.1.2. Varying Primary Constraints (ro e , T e , W+ /H+) 

The parameter space becomes 8-dimensional if variation 
in the primary constraints (n c , T e , and W + /H + ) is allowed. 
In this case we are interested in how composition is affected 
by allowing these primary constraints to take on values re- 
flecting a wide range of observations (Section 2.2). Such an 
exercise illustrates how the choice of constraints affects the 
baseline solution. As with the baseline solution already dis- 
cussed, we set / H + = 1 in all cases since the results are only 
mildly sensitive to / H + (Section 5.1.1). 

To this end, the grid was searched to find parameter com- 
binations consistent with the following six constraint com- 
binations: 
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= 60 cm 3 


T c 


= 2eV 
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= 30 



The numbers in bold indicate the values used thus far. In 
each case, only one of the three constraints is different from 
the original set (n c = 60 cm" 3 , T c = 2 eV, W+/H+ = 12). 
This modest sampling allows us to illustrate the dependence 
of plasma composition to each primary constraint. 

Figure 7 shows the best fit for each case. The parame- 
ter combination giving model results in best agreement with 
each set of the constraints is shown in the left panel of each 
bar chart. In each chart, the original baseline case has been 
normalized to one, and the other two cases are given relative 
to this value (see Table 1 to determine the actual value for 
each parameter). We have also generated sensitivity con- 
tours similar to those presented in Figures 4 and 5 for each 
of these new constraint combinations (not shown). From 
the /diff = 0.05 contours, we have derived the ranges over- 
plotted on the parameters in Figure 7. No such ranges exist 
in the cases of T c = 1 eV and T c = 3 eV because the best 
fits have a total fractional difference of 0.43 and 0.19, re- 
spectively. 

Water-group ion composition ratios are plotted in the 
right panel of each bar plot. Densities and temperatures 
have not been shown because we are primarily interested 
in the effect on composition. The top bar chart shows 
that neutral source rate, hot-electron temperature, and hot- 
electron fraction respond monotonically to n e . In the middle 
chart, the neutral source and the hot-electron temperature 
decrease, while the hot-electron fraction and transport time 
increase with T e . 

Electron density imposes the weakest change in composi- 
tion, and the W + /H + ratio is the strongest driver. In partic- 
ular, H30 + is the most abundant water-group ion only when 
W + /H + = 30. This is an important result because obtaining 



the W + /H + ratio from the CAPS data is more straightfor- 
ward than obtaining the separate, individual abundances 
of each water-group species [Wilson et al, 2008]. Increas- 
ing the W + /H + ratio requires a lower hot-electron temper- 
ature. One reason for this is that the impact ionization rate 
of H depends less on the hot-electron temperature than do 
the water-group ionization rates (Figure 8). Therefore, as 
the hot-electron temperature lowers, the W + /H + ratio in- 
creases. H30 + on the other hand, is the only water-group 
ion that does not require hot-electron impact ionization to 
thrive. In fact, it serves as a sink for both OH + and H20 + 
ions via charge exchange (Table 3). Thus, as required by 
the fit in Figure 7, the hot-electron temp drops to raise the 
W + /H + ratio and in the process increases the HsO + /W + 
ratio. 

6. Hot-electron (/ e h) Modulation 

Steffi et al. [2008] found that the interaction of two hot- 
electron populations orbiting Jupiter with System III and 
System IV periods are required to recover the Cassini UVIS 
temporal observations of the Io torus composition, both in 
terms of amplitude and rotational period. Delamere and 
Bagenal [2008] propose that the observed azimuthal electron 
density modulation at Saturn [Gurnett et al, 2007] is also 
caused by an azimuthally varying hot-electron abundance. 
Motivated by their research, we present compositional sensi- 
tivity to a magnetic-longitude-dependent hot-electron frac- 
tion. Figure 9 shows the response of the Enceladus torus to 
a prescribed sinusoidal hot-electron fraction: 

/eh(A mag ) = /e b h sclinc + 0.0016 sin(A mag ), (10) 

where f^- Behne = 0.0046 and A mag is magnetic longitude with 
arbitrary phase. The modulation amplitude is chosen so 
that n c modulates by roughly a factor of two (ss 40-80 cm" 3 ) 
according to the Saturnian kilometric radio emission analy- 
sis by Gurnett et al. [2007] (see their Figure 2). All other 
parameters are fixed at the baseline values listed in Table 1. 

The top panel of Figure 9 illustrates that T c and n e are 
in phase with /eh- Impact dissociation and ionization by 
hot electrons (Table 3) drive the W + /H + ratio out of phase 
with /eh- All quantities are plotted on the same logarithmic 
scale to show how linearly they respond to the hot-electron 
modulation. 

The middle panel shows the water-group compositional 
response to / e h (plotted on the same logarithmic scale 
as the top panel). The baseline composition ratios at 
Amag = 0° (and 360°) are given in Table 1. The 
two hemispheres 0°-180° and 180°-360° are defined by 
a higher- and lower-than-baseline / e h, respectively. As 
/eh increases in the 0°-180° hemisphere, the hierarchy be- 
comes OH+ > H 2 0+ > 0+ > H 3 0+. As /eh drops in the 
180° — 360° hemisphere, H20 + begins to dominate and + 
becomes the minor species (H 2 + > H 3 + > OH+ > + ). 

The bottom panel illustrates how strongly and how lin- 
early composition and thermal-electron parameters respond 
to /eh- The plot shows the normalized change in each pa- 
rameter with respect to variation in / e h- Points are missing 
at the middle and endpoints because the plotted function 
diverges there as A/ ch = f eh - / e b h aseline -> 0. The solid 
and dashed lines respectively indicate quantities in and out 
of phase with / e h- Quantities in phase with / e h (n e , T c , 
+ /W + , and OH + /W + ) have essentially fiat curves, indi- 
cating linearity with / e h- + /W + responds to perturbations 
3-4 times less strongly than the rest. Out-of-phase quanti- 
ties (W+/H+, H 2 0+/W+, and H 3 0+/W + ) also behave lin- 
early with the exception of H20 + /W , which always reacts 
weakly and practically not at all when f ch is low. HsO + /W + 
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and W + /H + are always driven more efficiently than one-to- 
one. 

Correlations and anti-correlations between electron den- 
sity and composition such as those presented in Figure 9 
have already been observed at Jupiter by Steffi et al. [2006] . 
In particular, they find that S + is in phase, while S 3+ is 180° 
out-of-phase with equatorial electron density modulation in 
the Io torus. 

7. Conclusions 

We have compared output from our model to the con- 
straints on thermal-electron density (n c ), thermal-electron 
temperature (T c ), and mixing ratio of water-group ions 
to protons (W + /H + ) by exploring the space spanned by 
the following four parameters: neutral source rate (N SIC ), 
hot-electron temperature (T h), hot-electron density (n e h = 
/ehtte), and radial transport time scale (r tra ns). 

Our important results are: 

1. For the constraint choices of n c = 60 cm~ 3 , T c — 2 
eV, and W + /H + = 12, we find the following limits on the 
parameters: 

1-5 < AW(10- 4 cm" 3 s- 1 ) < 2.7 

12 days < T trans 

0-3 < / eh (%) < 0.9 

100 < T ch /eV < 250. 
For a volume of 27r(4i?s)(2-Rs) 2 , the source rate can be 
scaled to give a mass source rate of 100 < N Brc /(kg s -1 ) < 
180. We find that / H + (the fraction of protons confined to 
the equator) is strongly coupled to the hot-electron popu- 
lation and has not been constrained by this study (Section 
5.1.1). The solution space can be compared with future 
measurements of the parameters {N BIC , / c h, Teh, T trms ) and 
composition mixing ratios. Upper limits on UV power em- 
anating from the Enceladus torus from neutral oxygen at 
1304, 1356, and 6300 A would be very useful. 

2. With the full water-based chemistry, photo- plus im- 
pact ionization is nearly as important as charge exchange at 
providing energy by way of fresh pickup ions (Figure 2). 

3. The water-based chemistry (particularly dissociative 
recombination) increases the neutral-to-ion ratio from 12 
(the Delamere et al. [2007] oxygen-based model) to ^ 40. 
Further, the Enceladus torus remains neutral-dominated 
even when the thermal-electron temperature approaches the 
temperature of electrons in Jupiter's Io torus (6 eV). 

4. The H30 + /W + ratio is directly correlated with 
the W + /H + ratio (Figure 7), implying that HaO + is 
strongly anti-correlated with H + (Section 5.1.2). This re- 
sult is important because obtaining the W + /H + ratio is 
more straightforward than obtaining individual water-group 
abundances from the CAPS data. However, Sittler et al. 
[2008] have shown, based on CAPS SOI data, that H 3 0+ 
dominates the water-group near the orbit of Enceladus (their 
Figure 15). The dominance of H3O 4 " seen in the CAPS data 
has not been obtained by our model with the given con- 
straints. It is likely (manuscript in preparation) that the lo- 
cal interaction of the corotating plasma with the Enceladus 
plumes may contribute significantly to the HsO + abundance 
[Cravens et al, 2009]. 

5. Hot electrons are necessary in the Enceladus to en- 
hance ionization torus but do not directly contribute more 
that 1% of the total electron population. 

6. Significant variation in composition can be driven by 
a small perturbation in the hot-electron population (Figure 
9). 

The sensitivity study presented here will be useful in a 
future interpretation of longitudinal and temporal observa- 
tions (e.g., Gurnett et al. [2007]) of the Enceladus torus, 
especially in the context of a modulating hot-electron den- 
sity. 



Acknowledgments. This study was supported by NASA 
award number NNX08AB17G under the Cassini Data Analy- 
sis Program. CHIANTI is a collaborative project involving the 
NRL (USA), RAL (UK), MSSL (UK), the Universities of Flo- 
rence (Italy) and Cambridge (UK), and George Mason Univer- 
sity (USA). Generous allocations of computing time provided 
by the OU Supcrcomputing Center for Education and Research 
(OSCER) at the University of Oklahoma permitted the high- 
resolution grid calculation. 



References 

Bagcnal, F., J. D. Sullivan, and G. L. Siscoe (1980), Spatial dis- 
tribution of plasma in the Io torus, Geophys. Res. Lett., , 7, 
41-44. 

Barbosa, D. D. (1994), Neutral cloud theory of the Jovian nebula: 

Anomalous ionization effect of superthermal electrons, Astro- 

phys. J., , 430, 376-386, doi:10. 1086/174413. 
Barbosa, D. D., F. V. Coroniti, and A. Eviatar (1983), Coulomb 

thermal properties and stability of the Io plasma torus, Astro- 

phys. J., , 274, 429-442, doi:10. 1086/161459. 
Bridge, H. S., et al. (1981), Plasma observations near Saturn 

- Initial results from Voyager 1, Science, 212, 217-224, doi: 
10.1126/science.212.4491.217. 

Bridge, H. S., et al. (1982), Plasma observations near Saturn 

- Initial results from Voyager 2, Science, 215, 563-570, doi: 
10.1126/science.215.4532.563. 

Burger, M. H., E. C. Sittler, R. E. Johnson, H. T. Smith, 
O. J. Tucker, and V. I. Shematovich (2007), Under- 
standing the escape of water from Enceladus, Journal of 
Geophysical Research (Space Physics), 112, 6219 — h, doi: 
10.1029/2006JA012086. 

Cravens, T. E., R. L. McNutt, J. H. Waite, I. P. Robertson, J. G. 
Luhmann, W. Kasprzak, and W.-H. Ip (2009), Plume iono- 
sphere of Enceladus as seen by the Cassini ion and neutral 
mass spectrometer, Geophys. Res. Lett., , 36, 8106 — h, doi: 
10.1029/2009GL037811. 

Delamere, P. A., and F. Bagenal (2003), Modeling vari- 
ability of plasma conditions in the Io torus, Journal of 
Geophysical Research (Space Physics), 108, 1276 — h, doi: 
10. 1029 /2002 JA009706. 

Delamere, P. A., and F. Bagenal (2008), Longitudinal plasma 
density variations at Saturn caused by hot electrons, Geophys. 
Res. Lett, , 35, 3107-4- , doi:10.1029/2007GL031095. 

Delamere, P. A., F. Bagenal, V. Dols, and L. C. Ray (2007), 
Saturn's neutral torus versus Jupiter's plasma torus, Geophys. 
Res. Lett, , 34, 9105-+, doi:10.1029/2007GL029437. 

Dols, V., P. A. Delamere, and F. Bagenal (2008), A multispecies 
chemistry model of Io's local interaction with the Plasma 
Torus, Journal of Geophysical Research (Space Physics), 113, 
9208-+, doi:10.1029/2007JA012805. 

Dougherty, M. K., K. K. Khurana, F. M. Neubauer, C. T. Russell, 
J. Saur, J. S. Leisner, and M. E. Burton (2006), Discovery of a 
Dynamic Atmosphere at Enceladus from Cassini Magnetome- 
ter Observations, in 37th Annual Lunar and Planetary Sci- 
ence Conference, Lunar and Planetary Institute Conference 
Abstracts, vol. 37, edited by S. Mackwell and E. Stansbery, 
pp. 1 585-4-. 

Esposito, L. W., et al. (2005), Ultraviolet Imaging Spectroscopy 
Shows an Active Saturnian System, Science, 307, 1251-1255, 
doi:10.1126/science.H05606. 

Farmer, A. J. (2009), Saturn in hot water: Viscous evolution of 
the enceladus torus, Icarus, In Press, Corrected Proof, — , doi: 
DOI: 10.1016/j.icarus.2009.02.031. 

Frank, L. A., B. G. Burck, K. L. Ackerson, J. H. Wolfe, and J. D. 
Mihalov (1980), Plasmas in Saturn's magnctosphere, J. Geo- 
phys. Res., , 85, 5695-5708, doi:10.1029/JA085iAllp05695. 

Gurnett, D. A., A. M. Persoon, W. S. Kurth, J. B. Groene, T. F. 
Averkamp, M. K. Dougherty, and D. J. Southwood (2007), 
The Variable Rotation Period of the Inner Region of Saturn's 
Plasma Disk, Science, 316, 442-, doi:10.1126/science. 1138562. 



X - 10 



FLESHMAN ET AL.: ENCELADUS TORUS MODEL 



Gurnett, D. A., ct al. (2004), The Cassini Radio and Plasma 
Wave Investigation, Space Science Reviews, 114, 395—463, doi: 
10.1007/sll214-004-1434-0. 

Gurnett, D. A., et al. (2005), Radio and Plasma Wave Obser- 
vations at Saturn from Cassini's Approach and First Orbit, 
Science, 307, 1255-1259, doi:10.1126/science.H05356. 

Hansen, C. J., L. Esposito, A. I. F. Stewart, J. Colwell, A. Hen- 
drix, W. Pryor, D. Shemansky, and R. West (2006), Ence- 
ladus' Water Vapor Plume, Science, 311, 1422-1425, doi: 
10.1126/science.ll21254. 

Jurac, S., and J. D. Richardson (2005), A self-consistent model 
of plasma and neutrals at Saturn: Neutral cloud morphology, 
Journal of Geophysical Research (Space Physics), 110, 9220- 
+, doi:10.1029/2004JA010635. 

Jurac, S., M. A. McGrath, R. E. Johnson, J. D. Richardson, V. M. 
Vasyliunas, and A. Eviatar (2002), Saturn: Search for a miss- 
ing water source, Geophys. Res. Lett., , 29(24), 240,000-1. 

Lewis, G. R., N. Andre, C. S. Arridge, A. J. Coates, L. K. Gilbert, 
D. R. Lindcr, and A. M. Rymer (2008), Derivation of den- 
sity and temperature from the Cassini-Huygens CAPS elec- 
tron spectrometer, Planetary Space Science, 56, 901-912, doi: 
10.1016/j.pss.2007.12.017. 

Lichtenbcrg, G., N. Thomas, and T. Fouchet (2001), Detection of 
S(IV) 10.51 fim emission from the Io plasma torus, J. Geophys. 
Res., , 106, 29,899-29,910, doi:10.1029/2001JA900020. 

Moncuquet, M., A. Lecacheux, N. Meyer- Vernet, B. Cecconi, and 
W. S. Kurth (2005), Quasi thermal noise spectroscopy in the 
inner magnctosphere of Saturn with Cassini/RPWS: Electron 
temperatures and density, Geophys. Res. Lett., , 32, 20 — h, 
doi:10. 1029 /2005GL022508. 

Persoon, A. M., D. A. Gurnett, W. S. Kurth, G. B. Hospo- 
darsky, J. B. Groene, P. Canu, and M. K. Dougherty (2005), 
Equatorial electron density measurements in Saturn's inner 
magnctosphere, Geophys. Res. Lett., , 32, 23,105 — h, doi: 
10.1029/2005GL024294. 

Persoon, A. M., ct al. (2009), A diffusive equilibrium model 
for the plasma density in Saturn's magnctosphere, Journal 
of Geophysical Research (Space Physics), 114, 4211 — h, doi: 
10.1029/2008JA013912. 

Richardson, J. D., and E. C. Sittler, Jr. (1990), A plasma density 
model for Saturn based on Voyager observations, J. Geophys. 
Res., , 95, 12,019-12,031, doi:10.1029/JA095iA08pl2019. 

Richardson, J. D., A. Eviatar, and G. L. Siscoc (1986), Satel- 
lite tori at Saturn, J. Geophys. Res., , 91, 8749-8755, doi: 
10.1029/JA091iA08p08749. 

Richardson, J. D., A. Eviatar, M. A. McGrath, and V. M. Va- 
syliunas (1998), OH in Saturn's magnctosphere: Observations 
and implications, J. Geophys. Res., , 103, 20,245-20,256, doi: 
10.1029/98JE01127. 

Rymer, A. M., B. H. Mauk, T. W. Hill, C. Paranicas, 
D. G. Mitchell, A. J. Coates, and D. T. Young (2008), 
Electron circulation in Saturn's magnctosphere, Journal of 
Geophysical Research (Space Physics), 113, 1201 — h, doi: 
10.1029/2007JA012589. 

Schippers, P., et al. (2008), Multi-instrument analysis of elec- 
tron populations in Saturn's magnctosphere, Journal of 
Geophysical Research (Space Physics), 113, 7208 — h, doi: 
10.1029/2008JA013098. 

Schrcicr, R., A. Eviatar, and V. M. Vasyliunas (1998), A two- 
dimensional model of plasma transport and chemistry in the 
Jovian magnetosphere, J. Geophys. Res., , 103, 19,901—19,914, 
doi:10.1029/98JE00697. 

Shemansky, D. E. (1988), Energy branching in the Io plasma 
torus - The failure of neutral cloud theory, J. Geophys. Res., 
, 93, 1773-1784. 

Shemansky, D. E., P. Matheson, D. T. Hall, H.-Y. Hu, and T. M. 
Tripp (1993), Detection of the hydroxyl radical in the Saturn 
magnetosphere, Nature, 363, 329-331, doi:10.1038/363329a0. 

Sittler, E. C, et al. (2005), Preliminary results on Sat- 
urn's inner plasmasphcrc as observed by Cassini: Compar- 
ison with Voyager, Geophys. Res. Lett., , 32, 14 — h, doi: 
10.1029/2005GL022653. 

Sittler, E. C, ct al. (2006), Cassini observations of Saturn's inner 
plasmasphcrc: Saturn orbit insertion results, Planetary Space 
Science, 54, 1197-1210, doi:10.1016/j.pss.2006.05.038. 

Sittler, E. C, et al. (2007), Erratum to "Cassini observations of 
Saturn's inner plasmasphere: Saturn orbit insertion results". 
[Planetary and Space Science 54 (2006) 1197 1210], Planetary 
Space Science, 55, 2218-2220, doi:10.1016/j.pss. 2006.11.022. 



Sittler, E. C, et al. (2008), Ion and neutral sources and sinks 
within Saturn's inner magnctosphere: Cassini results, Plane- 
tary Space Science, 56, 3-18, doi:10.1016/j.pss.2007.06.006. 

Sittler, E. C, Jr., K. W. Ogilvie, and J. D. Scuddcr (1983), Sur- 
vey of low-energy plasma electrons in Saturn's magnetosphere 
- Voyagers 1 and 2, J. Geophys. Res., , 88, 8847-8870, doi: 
10.1029/JA088iAllp08847. 

Steffi, A. J., F. Bagcnal, and A. I. F. Stewart (2004a), Cassini 
UVIS observations of the Io plasma torus. II. Radial varia- 
tions, Icarus, 172, 91-103, doi:10.1016/j.icarus.2004.04.016. 

Steffi, A. J., A. I. F. Stewart, and F. Bagenal (2004b), Cassini 
UVIS observations of the Io plasma torus. I. Initial results, 
Icarus, 172, 78-90, doi:10.1016/j.icarus.2003.12.027. 

Steffi, A. J., P. A. Delamere, and F. Bagenal (2006), Cassini 
UVIS observations of the Io plasma torus. III. Observations 
of temporal and azimuthal variability, Icarus, 180, 124-140, 
doi:10.1016/j.icarus.2005.07.013. 

Steffi, A. J., P. A. Delamere, and F. Bagenal (2008), Cassini 
UVIS observations of the Io plasma torus. IV. Modeling tem- 
poral and azimuthal variability, Icarus, 194, 153-165, doi: 
10.1016/j.icarus.2007.09.019. 

Thomas, N., F. Bagenal, T. W. Hill, and J. K. Wilson (2004), 
The Io neutral clouds and plasma torus, pp. 561-591. 

Tokar, R. L., ct al. (2006), The Interaction of the Atmosphere 
of Enccladus with Saturn's Plasma, Science, 311, 1409-1412, 
doi:10.1126/science,1121061. 

Tokar, R. L., et al. (2008), Cassini detection of water-group pick- 
up ions in the Enceladus torus, Geophys. Res. Lett., , 35, 
14,202-+, doi:10.1029/2008GL034749. 

Trainor, J. H., F. B. McDonald, and A. W. Schardt 
(1980), Observations of energetic ions and electrons in 
Saturn's magnctosphere, Science, 207, 421—425, doi: 
10.1126/science.207.4429.421. 

Williams, J. D., P. A. Delamere, F. Bagcnal, D. Reisenfeld, and 
R. Flcshman (2008), Water Group Composition Near the Or- 
bit of Enceladus, AGU Fall Meeting Abstracts, pp. B1384+. 

Wilson, R. J., R. L. Tokar, M. G. Henderson, T. W. Hill, M. F. 
Thomsen, and D. H. Pontius (2008), Cassini plasma spec- 
trometer thermal ion measurements in Saturn's inner mag- 
netosphere, Journal of Geophysical Research (Space Physics), 
113, 12,218-+, doi:10.1029/2008JA013486. 

Wolfe, J. H., J. D. Mihalov, H. R. Collard, D. D. McKibbin, 
L. A. Frank, and D. S. Intriligator (1980), Preliminary re- 
sults on the plasma environment of Saturn from the Pioneer 
11 plasma analyzer experiment, Science, 207, 403—407, doi: 
10.1126/science.207.4429.403. 

Young, D. T., et al. (2005), Composition and Dynamics of Plasma 
in Saturn's Magnctosphere, Science, 307, 1262—1266, doi: 
10.1126/science.ll06151. 



Bobby Flcshman, Homer L. Dodge Department of Physics and 
Astronomy, University of Oklahoma, 440 W. Brooks, Norman, 
Oklahoma 73019, USA (fleshman@nhn.ou.edu) 



FLESHMAN ET AL.: ENCELADUS TORUS MODEL 



X - 11 



Table 1. Model constraints (see Section 2.2 for references) 
and output for the best fit baseline solution. The fit is defined 
as the combination of parameters T e h, / e hi /h+> r trans, and 
iVsrc that gives the smallest total fractional difference between 
the data and the model output. Notice that the model output 
for the constraints (in bold) agrees with the data to at least 
two significant figures. No reaction leads to H3O in our set of 
reactions (appendix). 



Data (Constraints) 



Baseline Fit 



n e /cm 3 : 60 
T e /eV : 2.0 
W+/H+ : f2 



Neutral 
Densities (cm -3 ) 

n H 
™H 2 
n Q 
"OH 
«-H 2 
™H 3 

Ion/Electron 
Densities (cm -3 ) 



Teh/cV 
fch 

/h+ 

Ttrans/days 

7V src /cm~ 3 s 



160 
0.46% 
1.0 
26 

2.0E-4 



Mixing Ratios 



720 
< 1 
700 
770 
190 



O+/W+ : 0.15 

OH+/W+ : 0.30 

H2O+/W+ : 0.37 

H 3 0+/W+ : 0.18 

O+/H+ : 1.8 

W+/H+ : 12 

Ion/Electron 
Temperatures (eV) 



n e 


60 


T e 


2.0 


"eh 


0.28 




160 


n H + 


4.6 


T H+ 


4.0 




< 1 


T H+ 


6.5 


n o+ 


8.4 




38 


n o++ 


0.078 


T o++ 


35 


n OH+ 


17 


T OH+ 


39 


n H 2 0+ 


20 


T H 2 0+ 


42 


n H 3 0+ 


9.8 


T H 3 0+ 


42 
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Table 2. Baseline lifetimes for each species by mechanism in 
descending order of frequency (Rate = 1/t). Electron-impact 
ionization and photoionization include processes that are both 
ionizing and dissociative. A listing of lifetimes by reaction can 
be found in the appendix. 





Mechanism 


Rate (s~ 


- 1 ) 


t (days) 




Mechanism 


Rate (s~ 


- 1 ) 


t (days) 


H 


Charge Exchange 


7.3x10" 


8 


1.6xl0 2 


H+ 


Charge Exchange 


6.1x10" 


-6 


1.9x10° 




Impact Ionization 


9.6x10" 


9 


1.2xl0 3 




Radial Transport 


4.4x10" 


-7 


2.6X10 1 




Photoionization 


8.0x10" 


10 


1.5xl0 4 




Recombination 


5.1x10" 


-9 


2.3xl0 3 


H 2 


Impact Dissociaton 


7.6x10" 


7 


1.5X10 1 


H+ 


Charge Exchange 


2.6x10" 


-6 


4.4x10° 




Charge Exchange 


5.2x10" 


8 


2.2xl0 2 




Dissociative Recomb. 1.2x10 


-6 


9.6x10° 




Photoionization 


6.9x10" 


10 


1.7xl0 4 




Radial Transport 


4.4x10" 


-7 


2.6X10 1 




Photodissociation 


4.9x10" 


10 


2.4xl0 4 












n 


Charge Exchange 


6.6x10" 


8 


i 7 v 1 

1.1 AlU 


o+ 

w 


Charge Exchange 


3.1x10" 


-6 


q 7 v 1 nO 

O. 1 X 1U 




Impact Ionization 


2.6x10" 


-8 


4.5xl0 2 




Radial Transport 


4.4x10" 


-7 


2.6X10 1 




Photoionization 


2.3x10" 


9 


5.0x10 s 




Impact Ionization 


7.5x10" 


-9 


1.5x10 s 














Recombination 


1.9x10" 


11 


6.1xl0 5 














Radial Transport 


4.4x10" 


-7 


Z.D X lu 














Charge Exchange 


3.7x10" 


-7 


3 2x 10 1 














Recombination 


1.1x10" 


-10 


l.OxlO 5 


OH 


Photodissociation 


5.5x10" 


8 


2.1x 10 2 


OH+ 


Charge Exchange 


1.1x10" 


fi 


l.lxlO 1 




Impact Ionization 


3.5x10" 


8 


3.3xl0 2 




Dissociative Recomb. 5.8x10 


-7 


2.0X10 1 




Impact Dissociation 2.7x10 


8 


4.3xl0 2 




Radial Transport 


4.4x10" 


-7 


2.6X10 1 




Charge Exchange 


1.6x10" 


8 


7.4xl0 2 














Photoionization 


3.7x10" 


9 


3.1xl0 3 














Impact Dissociation 4.9x10 


7 


2.4X10 1 


H 2 0+ 


Dissociative Recomb. 1.2x10 


-6 


9.9x10° 




Charge Exchange 


3.6x10" 


7 


3.2X10 1 




Radial Transport 


4.4x10" 


-7 


2.6X10 1 




Photodissociation 


1.3x10" 


7 


9.3X10 1 




Charge Exchange 


4.0x10" 


-7 


2.9X10 1 




Impact Ionization 


5.2x10" 


8 


2.2xl0 2 














Photoionization 


4.5x10" 


9 


2.6xl0 3 






















H 3 Q+ 


Dissociative Recomb. 


8.1x10" 


-7 


1.4X10 1 














Radial Transport 


4.4x10" 


-7 


2.6X10 1 
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Table 3. List of the most important reactions for the base- 
line case. The left column gives the relevant impact, dissocia- 
tive, and photolytic reactions, and the right column gives all 
relevant charge exchanges and recombinations. The full set of 
reactions are given in the electronic appendix, but steady-state 
densities and temperatures are all within 3% of the properly 
calculated values when only the above reactions are turned on. 



Dominant Reactions 



H + c h 




H+ + 2c 


H+ + H 




H + 


H+ 


O + c h 




0+ + 2e 


H+ + 


O 




H + 


0+ 


OH + e 




OH+ + 2c 


H+ + 


OH 




H + 


OH+ 


OH + e h 




OH+ + 2c 


H+ + 


H 2 




H + 


H 2 0+ 


H 2 + c h 




H 2 0+ + 2e 


o+ + 


H 




O + 


H+ 


H 2 + c h 




OH+ + H + 2e 


o+ + 


O 




O + 


0+ 


H 2 + c h 




H+ + OH + 2c 


o+ + 


OH 




O + 


OH+ 


OH + c 




O + H + e 


o+ + 


H 2 




O + 


H 2 0+ 


OH + c h 




O + H + e 


OH+ 


f OH 




O + 


H 2 0+ 


H 2 + c 




OH + H + e 


OH+ 


f H 2 




OH 


+ H 2 0+ 


H 2 + c h 




OH + H + e 


OH+ 


f H a O 




O + 


H3O+ 


O + 7 




0+ + e 


H 2 0+ 


+ H 2 




OH 


+ H3O+ 


OH + 7 




OH+ + c 


H 2 0+ 


+ H 2 




H 2 + H 2 0+ 


OH + 7 




O + H 


OH+ 


f e 




+ H 


H 2 + 7 




H + OH 


H 2 0+ 


+ e 




OH 


+ H 


H 2 + 7 




H 2 + O 


H s O+ 


+ e 




OH 


+ H 2 
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2.0 xl(T 4 cm^ 3 s" 1 



125 kg s- 1 




e, O, 0+, 0++, 
H, H+, H 2 0, H 2 0+, 
OH, OH+, H 3 0, H 3 0+ 
H2, 




Figure 1. Particle flow for the baseline case. We have 
assumed a torus volume of 27r(4i?s)(2-Rs) 2 to calculate 
the volumetric mass flow. The percentages given here 
are for particle number (not mass). Individual species 
contributions to the particle outflow can be found in the 
electronic appendix. 




Figure 2. Energy flow for the baseline case. We have 
assumed a torus volume of 27r(4i?s)(2-Rs) 2 to calculate 
the volumetric energy flow. Individual species contribu- 
tions to the energy outflow can be found in the electronic 
appendix. 
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Figure 3. Sensitivity plots of fdis for every parameter 
combination. In each case, the remaining three parame- 
ters are fixed at their baseline values (Table 1) to ascer- 
tain trends due solely to variation of a single parameter 
at a time. The intersection of dashed lines indicates the 
baseline solution, and the gray shading inside of fdis = 1 
is intended to guide the eye. (Contours are plotted loga- 
rithmically.) 
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Figure 4. Sensitivity plots between hot-electron tem- 
perature (T e h) and hot-electron fraction (/ e h) gener- 
ated from the grid search. The over-plotted box rep- 
resents the baseline solution, discussed in Section 4. The 
/diff = 0.05 curve (Panel 1) has been over-plotted on 
Panels 5-13 in gray. The three parameters in Panels 5- 
7 (rtrans, N src , and /h+) have taken on values yielding 
best agreement between model output and n e =60 cm~ 3 , 
T c = 2 eV, and W+/H+ = 12. Combinations of T ch and 
/eh within the gray contour are consistent with these con- 
straints. All panels are discussed in Section 5.1. 
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Figure 5. Sensitivity plots between radial transport 
time scale (r tr ans) and hot-electron fraction (/ e h) gener- 
ated from the grid search. The over-plotted box repre- 
sents the baseline solution, discussed in Section 4. The 
/diff = 0.05 curve (Panel 1) has been over-plotted on 
Panels 5-13 in gray. The three parameters in Panels 5- 
7 (iV src , rtrans, and / H + ) have taken on values yielding 
best agreement between model output and n e = 60 cm -3 , 
T e = 2 eV, and W + /H + = 12. Combinations of rtrans and 
/ e h within the gray contour are consistent with these con- 
straints. All panels are discussed in Section 5.1. 
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Figure 6. Solutions for various proton dilution factors 
(/ H + = 0.7,0.85, 1.0) found from the grid search in Sec- 
tion 5.1.1. The parameters in the left panel are normal- 
ized to the values listed in Table 1. In the right panel 
we show the model output for water-group composition. 
Proton production is strongly controlled by impact ion- 
ization of hydrogen by hot electrons, so / H + is coupled to 
T ch and / e h- This coupling diminishes the significance of 
/ H + and effectively reduces the parameter space to four 
dimensions. 
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Figure 7. Solutions for primary constraint values rep- 
resenting a wide range of observations found from the 
grid search in Section 5.1.2. In each case, n e = 60 cm~ 3 , 
T c = 2 eV, and W + /H + = 12 unless otherwise speci- 
fied. The parameter fits given in the left panel of each 
bar chart are normalized to the fit given in Table 1. The 
over-plotted ranges on the parameters are derived from 
the corresponding /dig = 0.05 sensitivity contours. No 
such ranges exist for the T c = 1 eV and T c = 3 eV so- 
lutions because the best fits have a fractional difference 
of 0.43 and 0.19, respectively. The corresponding model 
output for water-group composition is presented on the 
right. The choice of electron density (n c ) has a weak 
effect, while the ratio W + /H + ratio alters composition 
markedly. 
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Figure 8. Electron-impact reaction rates for the water 
group (O + OH + H2O) and hydrogen. The bottom panel 
emphasizes that water-group reaction rates fall faster 
than hydrogen reaction rates as T ch increases. 
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Figure 9. Composition variability due to hot-electron 
modulation. The hot-electron fraction, / e h, alone drives 
the variation since all other parameters are held at their 
baseline values. The quantities + /W + , OH + /W + , n e , 
and T c are in phase while W+/H+, H 2 + /W + , and 
H30 + /W + are out of phase with / e h. The bottom panel 
shows how strongly and how linearly each quantity re- 
sponds to perturbations in / e h; the solid lines represent 
quantities in phase with / c h, and the dashed lines repre- 
sent quantities out of phase with / C h. 



